Universal properties in ultracold ion-atom interactions 
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We present some of the universal properties in ion-atom interaction derived from a newly for- 
mulated quantum-defect theory for — 1/r 4 type of long-range interactions. For bound states, we 
present the universal bound spectrum, namely the equivalent of the Rydberg formula, for ion-atom 
systems. For scattering, we introduce the concept of universal resonance spectrum to give a system- 
atic understanding of many resonances present in ion-atom scattering. The theory further provides 
a method for an accurate spectroscopic determination of the atomic polarizability. It also suggests 
the existence of atom-like molecules, in which multiple atoms orbit around a heavy ion. 

PACS numbers: 34.10.+x,34.50.Cx,33.15.-e,03.65.Nk 



As experimental techniques for preparing and manip- 
ulating cold atomic ions improve [l|, there is a rapidly 
growing interest in ion-atom interactions at ultracold 
temperatures 0-0] ■ While it is intuitively clear that 
an ion can exert a much stronger influence on its en- 
vironment than either an atom or a molecule, the pre- 
cise nature of this influence, especially at the quantum 
level, is far from being understood. For example, an ion- 
atom system is one of the few types of two-body systems 
for which their loosely-bound states, namely the highly 
excited rovibrational states right below the dissociation 
limit, are yet to be observed experimentally or fully char- 
acterized theoretically p. 

This work presents part of a newly formulated 
quantum-defect theory (QDT) for ion-atom interactions. 
It gives precise explanation and characterization of the 
meaning and the consequence of the "strong influence" 
of an ion on a neighboring atom. On a more technical 
level, it is a version of the QDT for —1/r 4 type of long- 
range potential 0-0] that brings our understanding of 
ion-atom interactions to the same level as, and in certain 
areas exceeding, our current understanding of ultracold 
atom-atom interactions 10]. The theory is a result of 
combinin g co nceptual developments, built further upon 



Refs. llHl3l|. with improved mathematical understand- 



ing of the modified Mathieu functions [14| - |l6l | , especially 
for negative energies. The latter development, which al- 
lows for the efficient determination of all QDT functions 
for —1/r 4 potential, is achieved by solving the modified 
Mathieu equation using techniques we have previously 



17J and l/r J 



developed for the analytic solutions of 1 /r 6 
[lif potentials. 

We focus here on the universal spectrum for ion-atom 
interactions. For bound states, it is the equivalent of the 
Rydberg formula for —1/r 4 type of long-range potential, 
formulated in a way to take advantage of the angular mo- 
mentum insensitivity of the short-range parameter that is 
a characteristic of atom-atom and ion-atom interactions 
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13l Il9j |. It allows for the determination of the en- 
tire rovibrational spectrum in the threshold region from 
a single short-range parameter, such as the quantum de- 



fect. More importantly, the concept of universal spec- 
trum is generalized here to positive energies to include 
scattering resonance positions. It facilitates a system- 
atic understanding of many resonances that are present 
in ion-atom scattering 0, 0, IH , within the same frame- 
work, and using the same short-range parameter, as the 
bound spectrum. A number of more specific results will 
be extracted out of the universal spectrum, including the 
energy bins II2I ] that contain various loosely-bound 
states or scattering resonances. Furthermore, we show 
that there exists a critical energy, below which all scat- 
tering resonances have positive widths that correspond to 
time delay [2(| , and above which all scattering resonances 
have negative widths that correspond to time advance 
[ioj ]. Such concepts, together with the concept of univer- 
sal width function to be briefly discussed there, provide 
what we believe to be the first systematic understanding 
of the global structure of scattering. Conceptually, they 
are expected to be generally applicable to all quantum 
systems with —1/r™ type of long-range potential with 
n > 2. We briefly outline selective qualitative features of 
ion-atom interaction before getting to the mathematical 
formulation and the results of the theory. 

Like other quantum systems around the threshold, 
ion-atom interaction in the ultracold regime is deter- 
mined primarily by its long-range interaction. Assuming 
that the atom involved has no permanent quadrupolc, 
such as the case for alkali-metal atoms in their ground 
states, the long-range ion-atom interaction has the form 
of V(r) ~ — C^/r 4 where C4 is related, in atomic units, 
to the static dipole polarizability of the atom, ot\, by 
C4 = cci/2. The qualitative features of ion-atom in- 
teraction can be understood through the length scale, 
p4 = (2^iC'4,/h 2 ) 1 / 2 , and other related scale parameters 
that are associated with the long-range potential. TableU 
gives samples of such parameters for alkali-metal sys- 
tems. The length scale /?4 is a measure of the size of 
the last few loosely-bound states. The corresponding en- 
ergy scale, se — (h 2 /2/i)(l/ p^) 2 , sets the scale for their 
energy spacing. The related time scale st — H/se sets 
the scale for, e.g., the lifetime of a shape resonance, and 
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TABLE I: Sample scale parameters for ion-atom systems. Ca, 
in atomic units, is given by ai/2 where a 1 is the static dipole 
polarizability of the atom. The Pa = {^^lCa/Ti 2 ) 1 ^ 2 is the 
length scale. The sb = (h 2 /2u)(1/Pa) 2 is the corresponding 
energy scale. It is given both in units of microkelvin and in 
units of kHz. 



System 


C4 (a.u.) 


Pa (a.u.) 


s E /k B (fj,K) 


s E /h (kHz) 


87 Rb+- 6 Li 


82.06" 


1297 


9.148 


190.6 


87 Rb+- 23 Na 


81.30 6 


2321 


0.8841 


18.42 


6 Li+- 6 Li 


82.06° 


948.5 


32.01 


667.0 


23 Na+- 23 Na 


81.30 b 


1846 


2.212 


46.08 


40 K + - 40 K 


145. 1 6 


3251 


0.4101 


8.545 


133 Cs +_133 Cs 


199.9 6 


6959 


0.02691 


0.5608 



"Using a\ from Ref. 
'Using 01 from Ref. 



7r/3| sets the scale for low-energy scattering cross section. 

A large length scale (3a is what characterizes the 
"strong influence" of an ion on its environment. From 
Table U it is clear that ion-atom interaction has a much 
longer length scale than either the atom-atom interac- 
tion, with a corresponding length scale of (3§ ~ 100 
a.u., or the electron-atom interaction, with a correspond- 
ing length scale of /?4 ~ 10 a.u. More importantly, 
the /?4 for ion- atom is hundreds times greater than rrj, 
the distance below which the ion-atom interaction starts 
to deviate from the — Cj/r 4 behavior, and which can 
be estimated to be of the order of ro ~ 10-30 a.u.. 
The implication of this large scale separation is that an 
ion-atom system typically supports thousands of bound 



states, and, from Levinson theorem [23}, thousands of 
scattering resonances, with many of them around the 
threshold following the universal behavior as described 
by the angular-momentum-insensitive formulation of the 
quantum-defect theory (AQDT) [HEl. 

Universal bound spectrum: In AQDT [ll|, [l3j], the 
bound spectrum of a two-body system with — l/r n [n > 
2) type of long-range potential is given rigorously by the 
solutions of 



c(n) / 
Xi (e s 



(1) 



Here the Xi function is a universal function of a scaled 
energy, e s = e/se, that is uniquely determined by the 
exponent of the long-range interaction n and the angular 
momentum I. All the short-range physics is encapsulated 
in the parameter K c (e,l). It is a short-range K matrix 
that is a meromorphic function of both e and I [l3| . The 
K c parameter is further related to the quantum defect 
[i c , which is defined to have a range of < [i c < 1 , by 
K c (e, I) = tan[7r^ c (e, I) + nb/2] where b = l/(n - 2) (nf. 

For n — 4, we have obtained all QDT functions of 
Ref. [13j through an improved understanding of the mod- 
ified Mathieu equation 14-16j]. In particular, we have 



obtained the universal Xi 



^ (4) =tan(7nV2) T 



function for n = 4: 



1 + M e 2 , 



(2) 



where v is the characteristic exponent for — 1 /r 4 type of 
potential [14jl6l l24j , and 



1- 2v \e 



/2 / r[l-(i/ + n,)/2] \ 
S ' \T[l + (v + v )/2)}) 
r[l-{v-v )/2}\ (C e ,i{-v) 



MM = 



T[l + {v-v )/2]J \C eJ (+v) 
in which i>o = / + 1/2, 

oo 

with Q(v) being given by a continued fraction 



Qiy) = 



1 



,(3) 



(4) 



(5) 




Quantum Defect 



FIG. 1: (Color online) The universal spectrum for two-body 
quantum systems with — 1/r 4 type of interactions, including 
both the universal bound spectrum for e < 0, where the y-axis 
is defined as — (— e s ) 1//4 , and the universal resonance spectrum 
for e > 0, where the y-axis is defined as (e 3 ) 1//4 . For any 
two-body quantum system with — 1/r 4 type of long-range po- 
tential, the bound spectrum and the resonance spectrum are 
given by the cross points between the curves plotted and the 
curves representing u c (e,l). For systems such as ion-atom, 
for which Pa S> J*o, M c (e, I) is approximately an /-independent 
constant, allowing the determination of the entire rovibra- 
tional spectrum, and the resonance spectrum, from a single 
parameter, as illustrated by a vertical line in the figure. 

Equations Jl]) and (J2J) give the universal bound spec- 
trum for —1/r 4 type of long-range interaction. Figure [T] 
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TABLE II: Energy bins for the first few partial waves. The 
i-th bound state of angular momentum I, with i = 1 cor- 
responding to the least-bound state, is to be found within 
B_i < e s < -B_i+i for i > 1, and within B_i < e s < for 
z = 1. Shape resonances of angular momentum I can only 
exist between < e s < B c . The i-th diffraction resonance 
of angular momentum defined as a resonance with negative 
width, is to be found within Bi-i < e s < Bi for i > 1. A 
zeroth diffraction resonance may exist within B c < e 3 < Bo, 
depending on the quantum defect. 



B x 


I = 


I = 1 


I = 2 


I = 3 


I = 4 


B 4 


4782.3 


10266 


19073 


32206 


50780 


B, 


1734.9 


4480.1 


9371.7 


17218 


28940 


B 2 


442.02 


1553.1 


3886.6 


8057.0 


14792 


Bi 


52.766 


350.37 


1199.8 


3021.3 


6348.0 


Bo 





20.514 


174.75 


687.38 


1896.0 


B c 





4.8878 


35.505 


120.16 


292.38 


B-i 


-105.81 


-336.54 


-753.97 


-1407.4 


-2345.8 


B- a 


-1179.9 


-2432.9 


-4284.6 


-6835.1 


-10184 


B-3 


-5207.5 


-8840.6 


-13708 


-19961 


-27750 


B.i 


-15308 


-23247 


-33279 


-45607 


-60433 



gives one of its possible graphical representations, in com- 
bination with the universal resonance spectrum to be dis- 
cussed next. 

Universal scattering properties: In AQDT, the single- 
channel scattering properties are determined from [13J 



Ki = tan 5i 



Z c S n) K c 



Z f 



P) ( 



z 



; (n) 
fs 



Zgi n) K C 



(6) 

where Z^y (e s , I) are again universal QDT functions that 
are uniquely determined by the long-range exponent n 
and the angular momentum I [13| . 

We use here the definition of scattering resonance po- 
sitions as energies at which sin 2 (5; = 1, namely the en- 
ergies at which the scattering cross section reaches its 
unitarity limit 24]. Such locations can be determined 
as the roots of the denominator in Eq. ([6]). Defin- 
ing a generalized Xi function for positive energies as 



7 C (") I 



X'f '(ts) = Zf^/Zgl n \ the resonance positions can be 
formulated in a similar fashion as the bound spectrum, 
as the solutions of 



For n = 4, we have obtained 



Xi 



(4). , /0 , l-(-l)'M5tan[^(^-^)/2] 
~ n(W/ ) 1 + RP J tan[^ - VQ )/2] 



(7) 



(8) 



in which vq, and M £g i are the same as those defined 
earlier in relation to the bound spectrum. 

Equations ([7]) and © give the universal resonance 
spectrum for — 1/r 4 type of long-range interaction. The 
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FIG. 2: (Color online) Illustrating definitions of energy bins, 
and the critical energy B c , using the case d wave. The "Scaled 
Energy" is defined as in Figure [T] The y-axis is defined as 
Xi = [tan - (xf) — irb/2]/-K, in which tan - (xf) is taken to 
be within a range of it of [7rfc/2, tt + nb/2). Figure [T]is made 
up of multiple such curves for different I, except that they are 
turned around for a better visualization of the spectrum. 



e > portion of Fig. Q] gives one of its possible graphi- 
cal representations. Similar to a bound spectrum, which 
describes, over a set of discrete energies, the rise of a 
phase from zero to a finite value at the threshold, a res- 
onance spectrum describes its subsequent evolution back 
towards zero. For a scattering resonance located at e s i, 
which is one of the solutions of Eq. (J7J , the width of the 
resonance, again derived from Eq. ([6]), is given by 



dx 



in) 



de. 



(9) 



It shows that while the position of a scattering resonance 
depends on the short-range parameter, the width of the 
resonance, as a function of the resonance position, fol- 
lows a universal behavior that is uniquely determined by 
the long-range interaction. We call Eq. the universal 
width function. 

There arc a number of conclusions that can be drawn 
from the universal spectrum. We give here a brief sum- 
mary, with more complete discussions to be given else- 
where, (a) For ion-atom systems, namely systems with 
Pi ^ ?*o j the universal spectrum allows the determination 
of the entire rovibrational spectrum and the entire reso- 
nance spectrum in the threshold region, including states 
of different I, using a single parameter, assuming ct\ is 
known accurately. This single parameter can be deter- 
mined from the measurement of either a single bound 
state energy [11(, or a single resonance position, (b) 
From measurements of two or more bound state ener- 
gies and I or resonance positions, and using the method 
of Ref. Uj and its generalizations, the universal spec- 
trum can give not only an accurate prediction of the 
entire spectrum, but also an accurate, spectroscopic de- 
termination of the atomic dipole polarizability, with no 
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knowledge required of the details of the atomic interac- 
tion at short distances 
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The accuracy that can be 
expected with such a measurement is conservatively esti- 
mated to be 10~ 6 , far better than the best achieved us- 



ing atomic interferometry 25| . (c) The measurements of 
bound state energies and resonance positions are further 
facilitated by the energy bin concept Io|, 12 1 embedded 
in the universal spectrum. They give ranges of energies 
in which a certain state is to be found, as explained and 
tabulated in Table Inland illustrated in Fig. [2] The earlier 
version of the concept [10 , [I2I has been generalized here 
to include bins for resonance states, (d) The properties of 
the xi functions show that there exists a critical energy 

for each I, B c (l), defined by d%? /de 



B c (l), Xi 1S a piecewise monotonically decreasing func- 
tion of energy with dxf ^ /de s < 0, which, from Eq. ©, 
implies that all resonances occurring in this region have 
positive widths corresponding to time delays [20l j . Above 
B c (l), xi evolves into a piecewise monotonically in- 
creasing function of energy with dx\ /de s > 0, implying 
that all resonances above B c (l) have negative widths cor- 
responding to time advances [20j. Such resonances will 
be called diffraction resonances [2J] to distinguish them 
from shape resonances, which have positive widths. We 
note that the incorporation of diffraction resonances into 
the resonance spectrum is essential to complete the phys- 
ical picture of the evolution of the phase back towards 
zero, (e) The universal spectrum gives the maximum 
number of shape resonances that can exist for a particu- 
lar I. As can be observed from Figure [TJ this number is 
one for I = 1 through 4. Results for other I will be pre- 
sented elsewhere, which show, e.g., that the minimum I 
that can support two shape resonances is I = 7 for ion- 
atom interactions, (f) The results here provide explicit 
illustrations and confirmations of the general properties 
of — 1/r™ (n > 2) type of interactions derived earlier, 
such as the breakdown of the semiclassical approxima- 
tion around the threshold (see Ref. (Hf and the refer- 
ences therein), and the conclusion that the least-bound 
state for an ion-atom system is either an s state or a p 
state, depending on the quantum defect (l9j . (g) The 
large size of an ion-atom bound state, as compared to 
the range of atom-atom interaction, implies that a heavy 
ion is likely to be able to bound multiple light atoms, es- 
pecially in different vibrational states (orbits), in which 
atom-atom interactions are likely to be negligible. Such 
states, which can be expected to be much more stable 
than the states proposed by Cote et al. [3|, are concep- 
tually similar to the doubly or multiply excited states of 
an atom, except that the latter states have much shorter 
lifetimes due to strong electron-electron interactions. It 
will be interesting to find out how many atoms can be 
bound in such states, and how their properties depend on 
their configurations, and parameters such as the atom-ion 



B c (i) 
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mass ratio, atomic statistics, and atom-atom scattering 
characteristics. 

In conclusion, we have shown that in addition to the 
more obvious aspect of having a large scattering cross 
section, the "strong influence" of an ion on a neighbor- 
ing atom is reflected in the presence of a large number 
of bound states and scattering resonances around the 
threshold, a situation that would normally complicates 
its description. Fortunately, the very same "strong in- 
fluence", as characterized by a large /J4, simultaneously 
ensures that all such states follow universal properties 
determined solely by the long-range potential. Further 
studies of few-body and many-body systems involving 
ions will help to reveal their important roles in chemical 
reactions and in catalysis. 

I thank Jun Ye, Alex Dalgarno, Chris Greene, Mari- 
anna Safronova, and Hossein Sadeghpour for helpful dis- 
cussions at various stages of this work. It was supported 
by NSF. 
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